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MODELING  OF  EQUILIBRIUM  GAS  ADSORPTION  FOR 
MULTICOMPONENT  VAPOR  MIXTURES 
PART  II 

1  .  INTRODUCTION 

Single  vapor  equilibrium  adsorption  isotherms  are  often  used  to  assess 
the  relative  efficiencies  with  whicu  adsorbents,  such  as  activated  carbon, 
remove  specific  vapors  from  air  in  air  purification  schemes.  Single  vapor 
adsorption  isotherms  can  often  be  predicted  at  equilibrium,  subject  to  some 
limitations,  from  the  physical  properties  of  the  adsorbate  vapors  by  tech- 
niques  based  upon  the  Dubinin-Polanyi  concept  of  affinity  coefficient. 

Past  studies  have  provided  vast  amounts  of  adsorption  isotherm  data  for 
single  vapor  adsorption  on  various  adsorbents.  However,  the  practical 
conditions  under  which  adsorbents  are  employed  are  usually  quite  different 
from  the  ideal  laboratory  conditions  under  which  the  single  vapor  Isotherms 
are  determined.  For  example,  several  adsorbate  species  are  usually  pre¬ 
sent.  The  vapor  species  to  be  adsorbed  may  not  be  exposed  to  the  adsorbent 
under  conditions  where  equilibrium  can  be  readily  attained.  The  study  of 
multicomponent  kinetic  and  equilibrium  adsorption  on  an  adsorbent  is  very 
important;  however,  kinetic  and  equilibrium  adsorption  studies  are  still  at 
a  very  preliminary  stage  while  mixed  vapor  adsorption  studies  are  more 
complex  and  time  consuming.  This  problem  can  be  alleviated  to  some  extent 
if  mixed  adsorption  isotherms  can  be  predicted  from  single  vapor  adsorption 
data. 

Although  there  has  been  considerable  study  involving  the  thermodynamic 
properties  of  adsorbates  on  adsorbents,  relatively  few  studies  have  con¬ 
sidered  adsorption  kinetics.  Most  of  the  studies  to  date^"?  have  dealt  with  the 
kinetics  of  single  vapor  adsorption  behavior.  Many  of  these  studies  have 
made  use  of  Wheeler's  approach  to  the  kinetics  of  gas  adsorption  by  beds  of 
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adsorbent  granules.  This  method  has  yet  to  be  successfully  applied  to 
multicomponent  systems. 

The  development  of  experimental  methods  for  determining  binary  adsorp¬ 
tion  isotherms  and  the  kinetics  of  binary  mixture  adsorption  was  completed 
in  the  first  year  of  the  project.  In  the  second  year,  the  adsorption  of 
CHCl^/CCljj,  n-Hexane/Benzene,  and  Cf^C^/CHCl^  binary  systems  on  BPL- 
activated  carbon  was  investigated.  BPL  is  a  designation  assigned  by  Calgon 
Corporation.  In  the  past  year,  additional  equilibrium  and  kinetic 
adsorption  studies  were  carried  out  on  n-hexane/benzene,  CH-^C^/CHCl^  and 
n-hexane/C^C^  binary  mixtures.  Also,  potentially  useful  theoretical 
models  were  applied  to  all  the  binary  systems  studied  during  the  reporting 
period.  Results  and  conclusions  from  these  studies  are  presented  in 
addition  to  the  recommendations  for  work  to  be  carried  out  in  the  future. 


BACKGROUND 


Equilibrium  Adsorption  of  Mixtures. 


Currently,  there  are  four  potentially  useful  theoretical  methods 
available  which  have  had  limited  success  in  predicting  the  adsorption 
characteristics  of  mixtures: 

a.  Dubinln-Polanyi  Pore  Filling  Theory1 

b.  John's  mixture  isotherm  model8 

c.  Myer's  ideal  adsorption  solution  theory^ 

d.  Proportionality  Method1® 


2.1.1 


Dubinin-Polanyi  Pore  Filling  Theory. 


The  Dubinln-Polanyi  theory  has  not  been  used  much  to  predict  multi- 
component  adsorption.  The  results  that  are  available  indicate,  however, 
that  the  theory  has  some  potential  for  application  to  multicomponent 
adsorption.  Bering  and  co-workers11*12  extended  the  Dubinin-Polanyi  equa¬ 
tion  to  tie  adsorption  of  mixtures  by  using  the  following  equation: 
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where  -  the  partial  molar  volume  of  mixture  component  i 

8^  -  the  partial  molar  affinity  coefficient  of  mixture  component  i 

pai  -  the  saturated  vapor  pressure  of  mixture  component  i 

Pi  -  the  equilibrium  pressure  of  mixture  component  i 

a^^  -  the  number  of  g  mole  of  component  i  adsorbed  per  gm  of  adsorbent 

-  the  mole  fraction  of  component  i  in  the  adsorbed  phase 

WQ  and  B  -  constants  character izing  the  adsorbent 

In  practice,  the  quantity  ZN^v^  can  be  found  from  the  phase  diagram  of 
the  volume  solution  assuming  a  liquid-like  adsorbate,  while  EN^  can  be 
found  simply  according  to  an  additive  scheme. ^ 

In  the  theory  of  micropore  filling  in  the  case  of  an  individual 
component,  a  normal  liquid  at  the  given  temperature,  existing  in 
equilibrium  with  its  saturated  vapor  at  the  pressure,  PQ,  is  selected  as 
the  standard  state.  In  the  case  of  multicomponent  adsorption,  it  is 
unclear,  a  priori,  whether  the  state  of  a  solution  whose  composition  is 
equal  to  the  composition  of  the  adsorbed  phase  or  the  state  of  a  solution 
existing  in  equilibrium  with  vapor  whose  composition  is  equal  to  the 
composition  of  the  equilibrium  vapor  above  the  adsorbed  phase  should  be 
selected  as  the  standard  state.  However,  studies1^  have  shown  that 
Equation  (t)  is  fulfilled  well  in  both  methods  of  selecting  the  standard 
state.  Selecting  such  standard  states,  we  can  rewrite  Equation  (1)  for  a 
binary  mixture  of  vapors  in  the  following  form: 

_OT^  1  7 

W  "  a12v12  "  Wo  (1°S  <£»  1  (2) 

s12 
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where 


h 


Ipi 


Zp 


si 


Equation  (2)  was  found  to  be  applicable  to  several  systems.1**  Other  methods 
based  on  the  Dubinin-Polanyi  approach  have  also  been  described  in  litera¬ 
ture  . 1 5 . 1 6 

2.1.2  John's  Mixture  Isotherm  Model. 

This  model,  which  was  developed  by  John  and  others,®  assumes  that  the 
single  vapor  isotherms  for  species  i  can  be  represented  by  the  following 
equation: 

loglog  P°  -  Cj  ♦  D,.  log  W°  (3) 

where  -  a  constant 
D.  -  a  constant 

1  i 

W°  -  amount  of  adsorbate  in  cm®/g  at  pressure,  pt 
pi  -  (Pi/P3l)  1C)N 

and  superscript  'o'  denotes  pure  component.  N  is  an  integer  between  2 
and  6. 

A  similar  equation  describes  the  binary  vapor  (components  1  and  2) 
adsorption  isotherm: 

log  log  ? ^2  m  +  ^12  W.j 2  ( 4 ) 

where  C12  -  ^1C1  ♦  Y2C2 

°12  “  Y1D1  +  y2D2 

W12  -  the  amount  of  mixed  adsorbate 
P12  “  ^Pl  +  P2^^Ps1  +  PsS^1®^ 

and  Y1  and  Y2  are  the  mole  fractions  of  components  1  and  2  in  the  gas 
phase. 
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The  constants  C1 ,  C2,  D1 ,  and  D2  can  be  obtained  from  the  single  vapor 
isotherms  and  then  used  to  calculate  W12,  in  Equation  (4),  using  known  or 
assumed  P^2  values. 

From  the  model,  the  micropore  volume,  WQ,  can  be  computed  for  a  single 
vapor  component  i  as  follows: 

log  W0  -  (log  log  10N  -  Cl)/Di  (5) 

Similarly  the  micropore  volume  in  the  binary  mixture  case  is  given  by: 

log  WQ  -  (log  log  10N  -  C12)/D12  (6) 

John  and  others®  have  shown  that  this  method  can  be  applied  to  binary 
and  ternary  systems  to  compute  their  total  adsorption. 

2.1 .3  Ideal  Adsorbed  Solution  Theory. 

The  ideal  adsorbed  solution  theory  has  been  used  to  predict  mixed  gas 
adsorption  with  some  success  at  low  coverages/low  relative  pressure.^'17 
The  method  assumes  that  the  adsorbed  phase  forms  an  ideal  solution  and 
involves  determining  the  'spreading  pressure'  for  the  single  vapor  iso¬ 
therms.  The  calculation  is  made  as  follows: 

a.  Obtain  the  single  vapor  isotherms  for  pure  components  in  terras  of 
the  amount  absorbed,  (cc/g)  versus  equilibrium  pressure  (torr). 

b.  The  spreading  pressure  for  these  adsorbates  is  calculated  as 
follows: 


ifA  f  **  n  .  o 

RT  "  I  o  dPi 

n 


(7) 


where  ir  -  spreading  pressure 

A  -  specific  area  of  adsorbent 

n  -  total  number  of  moles  in  adsorbed  phase/gm  of  adsorbent 
o 

Pi  ■  equilibrium  vapor  pressure  of  pure  component 


c.  Calculate  the  vapor  pressure  of  pure  components  at  constant 


spreading  pressure. 

d.  Describe  the  amount  adsorbed  (n[*)  at  these  vapor  pressures 
(calculated  in  step  c)  from  step  a. 

e.  Calculate  the  adsorption  equilibria  for  both  components  at  a  desired 
total  pressure,  P,  using  the  following  two  equations: 


pyi  -  p°n1  (8) 

py2  ■  p2n2  (9) 

Adding  these  two  equations,  the  expression  for  an  ideal  liquid  solution  is 
obtained.  ,p  _  p° 

^l*oo  (10) 

P1  *  P2 


The  vapor  phase  composition  is  given  by 

Y. 


(ID 


N1  is  the  mole  fraction  of  component  1  in  the  adsorbed  phase  and  Y1  is  the 
mole  fraction  of  component  1  in  the  gas  phase. 

f.  The  total  amount  adsorbed  is  obtained  by: 


I.jfi  02) 

n  o 
ni 

g.  Finally,  the  amount  of  each  component  adsorbed  from  the  gas 
mixture  is  given  by:  ni  -  nNA. 

The  complete  isobaric  composition  diagram  for  any  mixture  is  obtained 
by  repeating  the  above  calculation  for  different  values  of  the  spreading 
pressure. 


2.1.4  Proportionality  Method. 


Recently  Jonas,  et  al.10  suggested  that  the  adsorption  behavior  of  a 
mixture  can  successfully  be  predicted  using  a  simple  technique  called  the 
proportionality  method.  In  this  technique  vapor  concentrations  are 
expressed  in  terms  of  their  mole  fractions.  The  equation  can  be  written  as 
follows : 

W12  -  Y1W1  ♦  Y2W2  (13) 

where  W12  -  amount  of  mixed  adsorbate  (g/g) 

W1  -  the  amount  adsorbed  of  component  1  from  the  single  vapor  (g/g) 

W2  -  the  amount  adsorbed  of  component  2  from  the  single  vapor  (g/g) 

Y1  -  mole  fraction  of  component  1  in  the  gas  phase  of  a  mixture 

Y2  -  mole  fraction  of  component  2  in  the  gas  phase  of  a  mixture 

Jonas,  et  al.10  have  applied  this  method  to  binary  systems  with  some 
degree  of  success. 

2.2  Kinetics  of  Adsorption. 

Studies  on  the  kinetics  of  adsorption  on  activated  carbon  have  been 
previously  reported5-7  for  the  case  of  single  vapor  adsorption  but  have  yet 
to  be  extended  to  multicomponent  mixtures.  However,  many  of  the  mathema¬ 
tical  equations  and  kinetic  processes  which  describe  these  phenomena  for 
single  vapors  should  also  be  applicable  to  multicomponent  mixtures.  Of 
particular  interest  is  the  approach  taken  by  Wheeler1®’1^  which  has  been 
used  successfully  for  single  vapor  adsorption  kinetics.5-7  Wheeler's  equa¬ 
tion,  which  is  based  on  the  principle  of  mass  conservation,  can  be  written 
as  follows: 

t„  “  (we/coQ)Cwb  “  PB  q  ln(C0/Cx)/kv]  (14) 

where  CQ  -  the  inlet  gas  concentration  in  g/cm^ 
kv  -  the  first  order  rate  constant  in  min-1 
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pg  *  the  bulk  density  of  packing  in  g/cra^ 

Wg  -  the  kinetic  saturation  capacity  in  g/g 
-  the  bed  weight  in  gm 
Q  -  the  volume  flow  rate  in  cm^/min 
tb  -  the  breakthrough  time  in  min 
Cx  -  the  exit  gas  concentration. 

From  a  plot  of  tb  versus  Wb,  the  saturation  capacity  (We)  and  the 
first  order  rate  constant  (ky)  can  be  obtained.  By  setting  tb  of  Equation 
(13)  equal  to  zero  and  solving  for  Wb  one  obtains 

Wb  -  pBQ/kv  in  (C0/Cx)  -  Wc  (15) 

where  W„  is  identified  as  the  critical  bed  weight,  or  that  weight  of  carbon 

v 

just  sufficient  to  reduce  CQ  to  Cx  under  the  test  conditions. 

3.  WORK  OBJECTIVES 

The  ultimate  objectives  of  the  study  are  to  develop  methods  for  pre¬ 
dicting  the  adsorptive  behavior  of  mixed  gas  systems  on  activated  carbon 
adsorbents  from  a  knowledge  of  the  adsorptive  properties  of  the  pure  compo¬ 
nents  and  to  determine  applicability  of  the  Wheeler  equation  to  multi- 
component  kinetic  adsorption. 

The  work  has  been  divided  into  the  following  four  phases:  (1)  deter¬ 
mination  of  equilibrium  adsorption  isotherms  on  BPL-actlvated  carbon  for 
various  binary  vapor  mixtures  at  25 #C,  (2)  prediction  of  binary  vapor 
adsorption  isotherms  on  BPL-activated  carbon  for  comparison  with  experi¬ 
mentally  determined  isotherms,  (3)  determination  of  breakthrough  parameters 
through  adsorbent  beds  for  single  vapors  and  binary  vapor  mixtures,  and  (4) 
testing  the  applicability  of  Wheeler's  equation  to  binary  systems. 
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Another  objective  has  been  to  investigate  the  effect  of  adsorbate 
polarity  on  the  predictive  techniques  and,  also,  the  effect  of  adsorbate 
polarity  on  the  mixed  vapor  adsorption  data  obtained  through  the  use  of 
Wheeler's  equation.  The  mixtures  investigated  experimentally  to  date 
contained  (i)  two  non-polar  components  (NP-NP),  (ii)  two  weakly  polar 
components  (WP-WP)  and  (iii)  a  non-polar  component  and  a  weakly  polar 
component  (NP-WP). 

4.  PROCEDURES  AND  RESULTS 

4.1  Binary  Equilibrium  Adsorption. 

The  apparatus  and  experimental  procedures  used  for  measuring  the 
single  and  binary  vapor  adsorption  isotherms  have  been  described  in  a 
previous  report  by  Reucroft  and  others.20 

The  activated  carbon  was  a  Pittsburgh-activated  carbon,  type  BPL,  IE- 
30  mesh,  having  an  internal  surface  area  of  about  1000  n^/g  and 
approximately  80t  of  the  internal  surface  area  associated  with  pores  less 
than  20  8  in  diameter. 

The  experimental  procedure  (Method  B)  used  in  the  equilibrium  studies 
of  mixed  vapors  has  been  previously  described.21 

4.2  Results  and  Discussion. 

Figure  1  shows  the  single  vapor  adsorption  isotherms  for  various 
vapors  on  BPL  activated  carbon  at  room  temperature  in  terms  of  amount 
adsorbed  (gm/gm  of  activated  carbon)  versus  equilibrium  vapor  pressure. 

Figures  2-5  show  the  experimental  single  and  mixed  vapor  adsorption 
isotherms  of  the  CCl^-CHClj  binary  system  and  comparison  of  experimental  and 
theoretical  adsorption  Isotherm  data.  The  theoretical  adsorption  capacities 
were  calculated  using  D-P  theory  and  John's  equation.  Agreement  between  the 
predicted  and  experimental  equilibrium  capacities  is  quite  good.  Therefore, 
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Figure  2.  Single  Vapor  and  Mixed  Vapor  Adsorption  Isotherms 
(CCl^-CHClj  System)  on  BPL-Activated  Carbon 
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Figure  3.  Experimental  Mixed  Vapor  Isotherms  (CCl^(20)-CHCl3(80)  Mixture) 
Compared  with  Theoretical  Isotherms  on  BPL-Activated  Carbon 
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Experimental  Mixed  Vapor  Isotherms  (CCl4<80)-CHCl3(20)  Mixture) 
Compared  with  Theoretical  Isotherms  on  BPL -Activated  Carbon 
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Figure  6.  Single  Vapor  and  Mixed  Vapor  Isotherms  (Cfi  Wu  System)  on 
BPL-Activated  Carbon 
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Figure  8.  Experimental  Mixed  Vapor  Isotherms  (C^Hg (50)-C^H^4(50)  Mixture) 
Compared  with  Theoretical  Isotherms  on  BPL- Activated  Carbon 
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either  D-P  theory  or  John's  equation  can  be  used  to  calculate  adsorption 
capacity  of  CHCl^/CCl^  mixtures.  Experimental  isotherms  for  the 
n-hexane/benzene  system  are  shown  in  Figure  6.  Experimental  and  theoretical 
isotherms  are  compared  in  Figures  7  to  9  for  this  system.  Similar  results 
are  shown  for  the  n-hexane /CH2C12  system  in  Figures  10  to  14,  and  for  the 
CHCl^/Ci^C^  system  in  Figures  15  to  1 8. 

In  the  case  of  the  n-hexane/benzene  system,  the  experimental  isotherm 
agrees  well  with  John's  isotherm  equation  for  two  mixtures.  Better 
agreement  is  obtained  with  the  DP  theory  in  the  case  of  the  CgHg(25)~ 
c6hi 4(75)  Mixture. 

The  mixture  isotherms  are  very  similar  to  the  n-hexane  isotherms  for 
all  the  mixtures  investigated  in  the  case  of  the  n-hexane/ CH2C12  system. 

The  experimental  isotherms  do  not  agree  well  with  either  of  the  two  models 
for  this  system. 

In  the  case  of  the  CHClg/Cf^C^  system,  both  theoretical  models  tend 
to  underestimate  the  adsorption  amount  over  the  pressure  range 
investigated. 

4.3  Kinetics  of  Adsorption. 

The  experimental  procedure  and  a  schematic  of  the  binary  vapor  test 
apparatus  used  for  determining  the  kinetics  of  adsorption  were  described  in 
an  earlier  report  by  Reucroft,  et  al.20  The  concentration  of  each  component 
was  determined  from  the  traced  area  of  the  gas  chromatograph  peak.  The 
ratio  of  the  exit  concentration,  Cx,  to  the  inlet  concentration,  CQ,  was 
plotted  as  a  function  of  time,  t.  This  is  the  time  when  concentration 
ratio  reaches  0.01  (CX/CQ  ■  0.01)  is  called  breakthrough  time,  tb.  A 
linear  regression  analysis  of  the  data  for  tb  (break through  time)  vs  Wb 
(carbon  bed  weight)  for  a  binary  mixture  yielded  Wg  values  for  individual 
components  from  Wheeler's  equation.  Total  adsorption  of  the  bed,  Wm  was 
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Figure  10.  Single  Vapor  and  Mixed  Vapor  Isotherms  System) 

on  BPL-Activated  Carbon 
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Figure  11.  Experimental  Mixed  Vapor  Isotherms  (C6H14(9u)-CH2C12(10)  Mixture) 
Compared  with  Theoretical  Isotherms  on  BPL-Activated  Carbon 
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Figure  12.  Experimental  Mixed  Vapor  Isotherms  (CgHi4(75)-CH2Cl2(25)  Mixture) 
Compared  with  Theoretical  Isotherms  on  BPL-Activated  Carbon 
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Figure  13.  Experimental  Mixed  Vapor  Isotherms  (CgHi , (50)-CH2Cl2(50)  Mixture) 
Compared  with  Theoretical  Isotherms  on  BPL-Activated  Carbon 
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Figure  14.  Experimental  Mixed  Vapor  Isotherms  CgH]^(25)-CH2Cl2(75)  Mixture) 
Compared  with  Theoretical  Isotherms  on  BPL-Activated  Carbon 
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Figure  15.  Single  Vapor  and  Mixed  Vapor  Isotherms  (CHCI3-CH2CI2  System) 
on  BPL- Activated  Carbon 
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Figure  16.  Experimental  Mixed  Vapor  Isotherms  (CHC13(80)-CH2C12(20)  Mixture) 
Compared  with  Theoretical  Isotherms  on  BPL-Activated  Carbon 
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Figure  17.  Experimental  Mixed  Vapor  Isotherms  (CHCl3(50)-CH2Cl2(50)  Mixture) 
Compared  with  Theoretical  Isotherms  On  BPL-Activated  Carbon 
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Figure  18.  Experimental  Mixed  Vapor  Isotherms  (CHCl3(20)-CH2Cl2(80)  Mixture) 
Compared  with  Theoretical  Isotherms  on  BPL-Activated  Carbon 


also  obtained  by  directly  weighing  the  carbon  bed  before  and  after  the 
adsorption  process. 

4.4  Results  and  Discussion. 

A  typical  breakthrough  curve  for  all  the  single  vapors  shows  a 
sigmoidal  shape.  However,  some  of  the  components  in  the  mixtures  do  not 
follow  this  behavior.  All  the  investigated  binary  systems  show  that  the 
component  with  a  higher  saturated  vapor  pressure  is  displaced  regardless  of 
the  affinity  coefficient.  The  displaced  component  does  not  generally 
follow  a  sigmoidal  behavior.  The  displacement  becomes  evident  when  the 
Cx/C0  ratio  is  greater  than  1.0.  During  this  annual  reporting  period, 
additional  kinetics  studies  were  carried  out  for  n-hexane/benzene  (NP-NP) 
and  for  a  new  binary  system  CH2Cl2/n-hexane  (WP-NP).  In  addition,  four 
theoretical  models:  DP  theory,  John's  adsorption  isotherm  equation,  the 
ideal  adsorbed  solution  theory  and  the  proportionality  method  were  applied 
to  all  the  binary  mixtures  in  order  to  predict  the  adsorption  capacities. 
Predicted  values  were  then  compared  with  kinetic  and  gravimetric  capacities 
obtained  experimentally. 

Tables  1  and  2  show  linear  regression  analysis  of  the  tb  vs.  Wb  data 
for  the  two  binary  mixtures.  The  kinetic  capacities,  We,  obtained  from  the 
regression  analysis  and  listed  in  Tables  3  and  4  along  with  equilibrium 
capacities  obtained  gravimetrical ly  (Wg)  and  the  Wm  values  obtained  by 
direct  weighing. 

In  the  case  of  the  n-hexane/benzene  mixture  agreement  between  the 
three  adsorption  capacities  (We,  Wg  and  Wm  )  values  is  quite  good. 

However,  in  the  case  of  the  CHgClg/n-hexane  (WP-NP)  binary  system,  the 
kinetic  experimental  values,  We,  is  considerably  higher  than  both  the 
capacity  obtained  by  direct  weighting  (Wm)  and  the  equilibrium  gravimetric 
capacity  (Wg).  In  the  two  binary  systems  (CHCl^/CClj,,  CHgC^/CHCl^) 


studied  last  year,  it  was  also  observed  that  We  values  were  always  signifi¬ 
cantly  higher  than  the  Wg  and  Wm  values.  It  is  clear  that  Wheeler's 


Table  1.  Breakthrough  Time  (tb)  as  a  Function  of  Bed  Weight  (Wb)  for  the 
Ci-^C^/n-Hexane  Binary  System  on  BPL-Acti vated  Carbon.  (Flow 
Rate  -  400  cm^/min.,  Pj  »  25  torr,  T  -  298°K). 


Mole 

Bed 

Break through 

Regression 

Correlation 

Fraction 

Weight 

Time 

tb  (min) 

Equation 

Coefficient 

ch2ci2 

wb(g) 

ch2ci2 

n-Hexane 

0.20 

0.7000 

2.98 

2.99 

n- Hexane: 

1  .000 

1 .0007 

4.00 

4.37 

tb  -  4.434  Wb  - 

0.101 

1.1999 

5.00 

5.20 

CH2C12: 

tb  -  3.989  Wb  - 

0.137 

0.994 

0.50 

0.7002 

4.3 

5.0 

n-Hexane: 

0.998 

0.8517 

5.1 

6.0 

tb  -  7.317  Wb  - 

0.159 

1 .0008 

6.1 

7.2 

CH2C12: 

tb  -  5.926  Wb  - 

0.073 

0.998 

0.80 

0.7003 

4.0 

9.4 

n-Hexane : 

1  .000 

0.8505 

4.98 

11.8 

tb  -  15.542  Wb 

-  1.475 

1 .0007 

6.5 

14.0 

CH2C12: 

tb  -  5.29  Wb  - 

1  .893 

0.997 

equation  is  not  quite  effective  when  it  is  applied  to  weakly  polar  compo¬ 
nents  in  a  mixture. 

The  mole  fractions  of  the  more  volatile  components  in  the  adsorbed 
phase  were  calculated  from  the  WQ  values.  These  values,  as  a  function  of 
the  mole  fraction  of  the  same  component  in  the  gas  phase  are  shown  in 
Figures  19  to  21  for  the  three  mixtures  n-hexane/benzene,  Ci^C^/CHCl^  and 
n-hexane/C^C^.  It  appears  that  the  adsorbed  phase  mole  fraction  of 
CH2C12  in  the  CH2C12 /n-hexane  binary  system  is  similar  in  value  to  the  mole 
fraction  of  CHgC^  in  the  gas  phase,  except  at  70%  or  higher  concentration. 
However,  in  the  case  of  the  n-hexane/benzene  binary  system,  the  mole  frac¬ 
tion  of  n-hexane  in  the  adsorbed  phase  is  similar  in  value  to  the  mole 
fraction  of  n-hexane  in  the  gas  phase  only  at  low  concentration  (<0.3).  In 
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Vapor-Adsorbed  Phase  Equilibria  for  n-Hexane /Benzene  Binary 
Mixture  (BPL-Activated  Carbon,  T  *  25°C,  Py  *  25  torr,  flow 
rate  *  400  cm^/min.) 


in  Vapor 


NcH2C,2m  AdsorbeC*  Phas® 


Figure  20.  Vapor- Adsorbed  Phase  Equilibria  for  CHCWCHoCl,  Binary 
Mixture  (BPL-Activated  Carbon,  T  •  25°C,  PT  «  25  torr, 
flow  rate  ■  400  cm3/rain.) 
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Figure  21.  Vapor-Adsorbed  Phase  Equilibria  for  n-Hexane/Cl^C^  Binary 
Mixture  (BPL -Activated  Carbon,  T  -  25°C,  PT  ■  25  torr, 
flow  rate  ■  A00  cm^/min.) 
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Table  2.  Breakthrough  Time,  tb,  as  a  Function  of  Composition  and  Bed 
Weight  (Wb)  for  n-Hexane/Benzene  Binary  Mixtures* 


Mole  Fraction4, 

th  (min) 

wb  (g) 

Regression 

n-Hexane 

n-Hexane 

Benzene 

0.0 

- 

5.15 

0.8000 

Benzene : 

- 

6.00 

0.9037 

tb  -  7.7743Wh-1 .0403 

— 

7.55 

1 .1006 

9.85 

1  .4033 

corr.  -  0.9999 

0.2 

4.50 

4.50 

0.8004 

Benzene : 

5.95 

5.95 

1 .0020 

tb  -  7.4478Wb-.0839 

7.40 

7.50 

1  .2032 

corr.  -  .998 

n-Hexane:  U-7.1  996Wb~-1  .2630 
corr.  »  1.0000 

0.35 

3.80 

4.20 

0.7002 

Benzene : 

4.60 

4.85 

0.8509 

tb  -  6.27Wb-0.3589 
corr.  -  0.997 

5.50 

5.80 

1  .0003 

7.55 

8.50 

1 .4000 

n-Hexane:  tfc  -  5.37Wb-0.06l0 
corr.  •  0.9996 

0.50 

3.75 

4.10 

0.7005 

Benzene: 

4.45 

4.50 

0.8515 

tb  -  6.l8l3Wb  -  0.6835 

5.55 

5.85 

1 .0018 

corr.  ■  0.998 

6.65 

6.90 

1 .2001 

n-Hexane : 

8.10 

8.40 

1 .4097 

tb  -  6.2459Wb  -  0.5003 

0.80 

3.70 

3.20 

0.7019 

Benzene : 

4.45 

4.15 

0.8502 

tb  -  1 0 . 1 1 61  Wb  -  4.1374 

7.10 

9.20 

1.3109 

corr.  -  0.9961 

• 

n-Hexane : 

tb  -  5.6232Wb  -  0.2830 

corr.  -  0.9997 

1  .00 

4.2 

- 

0.5997 

n-Hexane: 

4.5 

- 

0.6875 

tb  -  8.l474Wb  -  1.049 

5.9 

- 

0.8498 

corr.  ■  0.9998 

7.1 

- 

0.9934 

10.35 

1.4028 

*BPL-Activated 
4Vapor  Phase 

Carbon,  Flow 

Rate  - 

400  cm^/min, 

25 °C,  and  *  25  torn. 

Table  3-  Kinetic  Saturation  Capacities  (Wg),  Direct  Weight  Saturation- 
Capacities  (Wm)  and  Equilibrium  Capacities  (Wg)  for  the 
CH2Cl9/n~Hexane  Binary  System  on  BPL-Activatea  Carbon.  (Flow 
rate  -  400  cm^/min.,  Py  -  25  torr,  T  »  298°K) 


Mole 

Fraction 

W 

e 

W5 

g/g 

Wg 

wm 

m 

ch2ci2 

ch2ci2 

c6H14 

Total 

g/g 

g/g 

0.2 

0.0364 

0.1645 

0.2009 

0.270 

0.298 

0.5 

0.1370 

0.1690 

0.3060 

0.2725 

0.252 

0.8 

0.3030 

0.1442 

0.4472 

0.290 

0.266 

Table  4.  Kinetic  Saturation  Capacities  (W  ),  Direct  Weight  Saturation 

Capacities  (Wm)  and  Equilibrium  Capacities  (W_)  for  the  CgHg/n- 
Hexane  Binary  System  on  BPL-Activated  Carbon  TFlow  rate  -  400 
cm3/min.,  PT  -  25  torr,  T  -  298°K) 


Mole 

We 

we 

Wg 

wm 

Fraction 

g/g 

n-hexane 

g/g 

Total 

g/g 

g/g 

c6H14 

c6h6 

0.0 

— 

0.330 

0.330 

0.370 

0.328 

0.2 

0.067 

0.250 

0.318 

0.318 

0.335 

0.5 

0.130 

0.144 

0.275 

0.312 

0.312 

0.8 

0.208 

0.085 

0.296 

0.302 

0.285 

1 .0 

0.226 

— 

0.226 

0.270 

0.252 

the  case  of  the  CHgCl 

2/chci3 

binary  system 

(Nca  ci  - 

0.35,  0.5 

and  0.65 

data  were  reported  in 

the  last  year’s  report)  the  mole 

fraction 

of  ch2ci2 

in  the  adsorbed  phase 

is  less 

than  the  gas  phase  value 

for  the  whole  range 

of  concentration  studied.  It  appears  that  each  system  behaves  differently. 
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As  mentioned  previously,  all  four  theories  have  been  applied  to  the 


CHCl^/CCljj,  n-hexane/benzene,  Ci^C^/CHCl^  and  n-hexane/CH2Cl2  binary 
mixtures  in  order  to  predict  the  adsorption  capacities.  The  results  are 
shown  in  Tables  5*8  for  all  four  models  investigated.  These  tables  also 
compare  the  kinetic  adsorption  capacities  calculated  from  Wheeler's 
equation.  The  maximum  deviation  of  the  predicted  capacities,  with  respect 
to  Wg,  for  all  Binary  mixtures  except  n-hexane/CH2Cl2,  ranges  from  -7.0?  to 
19?.  In  the  case  of  n-hexane/CH2Cl2,  the  comparison  indicates  that  none  of 
the  models  are  particularly  successful  in  predicting  the  correct  amount 
adsorbed  (the  deviation  ranges  from  19?  to  91?).  Also  it  is  not  possible 
to  employ  the  ideal  adsorbed  solution  theory  for  this  particular  binary 
at  the  given  experimental  conditions.  It  appears  that  a  large  saturated  vapor 


Table  5.  Comparison  of  Adsorption  Capacities  Obtained  from  Several  Predictive  Models 
with  the  Experimental  Values  (We  &  W  )  for  the  Benzene/n-Hexane  Mixture 
(PT  -  25  torr,  flow  rate  -  400  cm^/mln) 


Mole 

Predicted  Sorption  g/g 

Fraction 

DP-Model 

John's  Model 

Myers '  Model 

Proportionality 

Method* 

We 

Wg 

n-Hexane 

Wmix  *  Dev+ 

wmix  *  Dev+ 

wmix  *  Dev+ 

Wmix  ?  Dev 

g/g 

g/g 

0.0 

0.330 

0.370 

0.20 

0.340  -7.0? 

0.326  -2.5? 

0.337  -5.3? 

0.319  0.3? 

0.318 

0.318 

0.50 

0.318  -1.9? 

0.314  -0.6? 

0.323  -3.5? 

0.307  1.6? 

0.275 

0.312 

0.80 

0.295  2.3? 

0.300  -0.6? 

0.312  3.3? 

0.296  2.0? 

0.295 

0.302 

1 .00 

0.226 

0.270 

*wmlx  “  ^Benzene  x  wBenzene^  +  ^Nn-Hexane  x  Wn-Hexane^ 
♦With  respect  to  Wg 


Table  6.  Comparison  of  Adsorption  Capacities  .Obtained  from  Several  Predictive  Models 
with  the  Experimental  Values  (We  &  Wg)  for  the  CHCl^/CCl^  Mixture  (Pj  • 

25  torr,  flow  rate  -  400  cm^/min) 


Mole 

Predicted  Sorption  g/g 

Fraction 

DP-Model 

John's  Model 

Myers'  Model 

Proportionality 

Method* 

we 

Wg 

CHCI3 

Wmix  *  Dev+ 

wmix 

%  Dev+ 

wmix  *  Dev+ 

wmix 

%  Dev+ 

g/g 

g/g 

0.0 

0.696 

0.669 

0.20 

0.636  1.9% 

0.609 

-2.4% 

0.658  5.4% 

0.616 

-1 .2% 

0.786 

0.624 

0.50 

0.621  2.1% 

0.590 

-3.0% 

0.636  4.6% 

0.599 

-1 ,5% 

0.796 

0.608 

0.80 

0.607  4.3% 

0.569 

-2.2% 

0.605  4.0% 

0.582 

0 

0.722 

0.582 

1.00 

0.781 

0.585 

*Wmix  *  (NCHC1  x  WCHC1  ^  +  (NCC1.  x  WCC1.  > 
3  3  Hi* 


♦With  respect  to  Wg 

pressure  (PQ)  difference  and  large  difference  in  affinity  coefficient  (B) 
between  the  two  components  plays  an  important  role  in  determining  the 
binary  vapor  adsorption  characteristics.  Table  9  lists  the  difference  in 
affinity  coefficient  and  saturated  vapor  pressure  for  the  binary  systems. 

In  addition  to  the  total  amount  of  a  mixture  adsorbed  (W12),  it  is 
also  important  to  be  able  to  predict  the  amount  of  each  binary  component 
that  contributes  to  W12.  A  semi-empirical  formula  suggested  by  Lewis'^  can  be 
used,  in  principle,  to  predict  the  amount  of  each  binary  component 
adsorbed: 

W,  /W|  +  W2/W|  -  1 

where  Wj  ♦  W2  -  W12  and  Wi ,  W2,  W12  are  the  amounts  adsorbed  of  component 
1,  component  2  and  the  mixture  at  pressures: 


P|  (note  P12 


12 


?1° 


Pi  *  P?). 


There  is  also  the  proportionality  method  as  mentioned  previously  that  can  be 
used  to  predict  the  amount  of  each  binary  component  adsorbed: 

*12  "  *1*1  +  *2*2 


These  two  methods  were  applied  to  calculate  the  adsorption  capacities 
of  individual  components  at  a  total  pressure  of  25  torr  and  at  different 
compositions  for  the  CHgC^/CHC^,  n-hexane/CH2Cl2  and  the  n-hexane /benzene 
binary  systems.  Comparison  of  the  predicted  and  the  experimental 
adsorption  capacities  are  shown  in  Tables  10-12. 


Table  7.  Comparison  of  Adsorption  Capacities  Obtained  from  Several 

Predictive  Models  with  that  of  Experimental  Capacities  (Wg  and 
Wg)  for  the  CHCI3/CH2CI2  Mixture  (PT  -  25  torr,  flow  rate  - 
400  cm3/min) 


Predicted  Sorption,  g/g 


Mole 

Fraction  DP-Model  John's  Model  Myers'  Model  Proportionality 

Method* 


CH2C12 

*mix 

*Dev  + 

wmlx 

*Dev  + 

wmix 

*Dev  + 

wmix 

*Dev  + 

g/g 

6 

g/g 

0.0 

— 

— 

— 

— 

— • 

— 

— 

0.692 

0.585 

0.20 

0.549 

0.9* 

0.568 

-2.5* 

0.541 

2.4* 

0.562 

-1.4* 

0.933 

0.554 

0.50 

0.486 

8.8* 

0.521 

2.3* 

0.498 

6.6* 

0.515 

3.4* 

0.950 

0.533 

0.80 

0.405 

19.4* 

0.471 

6.4* 

0.453 

9.9* 

0.467 

7.2* 

0.722 

0.503 

1  .00 

— 

— 

— 

— 

— 

— 

— 

— 

0.478 

0.323 

*  wmix  "  (N  x  W  )  ♦  (N  ♦  W  ) 

CH2C12  CH2C12  CHCI3  CHCI3 

*  *  deviation  with  respect  to  W 

© 


Table  8.  Comparison  of  Adsorption  Capacities  Obtained  from  Several 

Predictive  Models  with  the  Experimental  Values  (We  and  Wg)  for 
the  CHpC^/n-Hexane  Mixture  (PT  -  25  torr,  flow  rate  - 
400  cnw /min) 


Predicted  Sorption,  g/g 

Mole 


Fraction  DP-Model 

John ' s 

Model 

Myers' 

Model 

Proportionality 

Method* 

Wg 

g/g 

ch2ci2 

"mix 

%Dev+ 

"mix 

?Dev+ 

"mix 

^Dev* 

"mix 

*Dev+  g/g 

0.2 

0.330 

22.9 

0.512 

90.0 

— 

— 

0.319 

19.1  0.201 

0.270 

0.5 

0.372 

27.8 

0.507 

74.2 

— 

— 

0.354 

21.5  0.306 

0.273 

0.8 

0.385 

32.8 

0.497 

71 .3 

_ 

__ 

0.388 

33.7  0.447 

0.290 

x  w  ^  +  (Nn-Hexane  x  Wn-Hexane 
ch2ci2  ch2ci2 


mix 

With  respect  to  Wg. 


) 


Table  9.  Affinity  Coefficient  Difference  and  Saturated  Vapor  Pressure 
Difference  of  Components  in  Mixtures  Investigated 


Mixture 

Affinity  Coeff 
Difference 

Saturate 

Vapor  Press  Diff. 

cci„-chci3  (NP-WP) 

0.11 

83  torr 

chci3-ch2ci2  (WP-WP) 

0.39 

205 

C6H6~c6H14  (NP-NP) 

0.28 

54.6 

Cgf^  ^-C^Clg  (NP-WP) 

0.70 

249.5 

c6h14-ch3coch3  (NP-SP) 

0.61 

175 

Table  10  shows  that  the  maximum  percentage  deviations  of  the  capacity 
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From  Lewis's  Equation 
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From  Lewi's  Equation 


From  Lewis's  Equation 


using  Methods  1  and  2  for  n-hexane  are  26?  and  30?  respectively.  However, 
the  capacities  calculated  for  benzene  in  this  mixture  using  the  same  models 
gave  small  percentage  deviation.  In  the  case  of  the  C^C^/CHCl^  binary 
system,  both  the  methods  showed  very  poor  agreement  with  the  experimental 
data.  For  the  (^C^/n-hexane  mixture,  none  of  the  methods  were  successful 
in  predicting  the  amount  adsorbed  of  each  individual  component.  Deviations 
range  from  4?  to  80%.  The  dichloromethane/n-hexane  mixture  has  the  largest 
difference  in  properties  compared  to  the  other  systems.  Preliminary 
studies  on  the  acetone/n-hexane  mixture  also  indicate  similar  behavior. 

For  these  particular  types  of  mixtures,  the  mixture  with  a  higher  affinity 
coefficient  (n-hexane)  seems  to  dominate  the  adsorption  behavior. 

From  the  studies  on  these  four  mixtures,  it  has  become  apparent  that 
each  mixture  behaves  differently.  NP-NP  mixtures  are  ideal  for  studies. 
Most  of  the  theoretical  models  are  successful  in  predicting  the  adsorption 
isotherms  for  NP-NP  mixtures.  However,  deviation  from  theory  becomes 
greater  as: 

1)  Polarity  is  introduced  into  the  mixture. 

2)  Difference  in  affinity  coefficient  of  the  two  components  is  large. 

3)  There  is  a  large  difference  in  saturated  vapor  pressure. 

5.  RECOMMENDATIONS  FOR  FUTURE  WORK 

Efforts  within  the  last  two  years  have  dealt  with  detailed  kinetic  and 
equilibrium  adsorption  studies  of  CHClj/CClj,  (WP-NP),  CHCl^/CHgC^  (WP-WP), 
CgHllt/CgHg  (NP-NP)  and  CHgC^/CgHjjj  (WP-NP)  binary  mixtures  and  single 
vapors  on  BPL-activated  carbon.  In  addition,  preliminary  results  have  been 
obtained  in  the  case  of  acetone/n-hexane  and  acetone-dichloromethane  binary 
mixtures  and  single  vapors. 


For  the  next  phase  of  the  project,  several  objectives  have  been  formu¬ 
lated. 

1)  Continuation  of  binary  vapor  adsorption  studies  on  acetone/hexane  and 
acetone/CH2Cl2  binary  mixtures. 

2)  Additional  experimental  studies  on  strongly  polar-strongly  polar  binary 
mixtures  adsorbed  on  BPL-activated  carbon. 

3)  Verification  of  the  kinetic  saturation  capacity  (Wg)  obtained  by  the 
Wheeler  equation.  Wg  should  be  independent  of  bed  weight.  This  will  be 
evaluated  by  conducting  studies  on  a  wider  range  of  bed  weight.  Experiments 
such  as  these  may  shed  light  on  the  observation  of  high  Wg  values  in  some 
systems. 
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